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We present the first measurement of the angle (j>3 of the Unitarity Triangle using a model- 
independent Dalitz plot analysis of B ± — > DK ± , D — } Kgir + iT~ decays. The method uses, as 
input, measurements of the strong phase of the D — > K%-k + -k^ amplitude from the CLEO col- 
laboration. The result is based on the full data set of 772 x 10 6 BB pairs collected by the Belle 
experiment at the T(4S) resonance. We obtain (f>3 = (77.31^'g ± 4.1 ± 4.3)° and the suppressed 
amplitude ratio r_e = 0.145 ± 0.030 ± 0.010 ± 0.011. Here the first error is statistical, the second 
is the experimental systematic uncertainty, and the third is the error due to the precision of the 
strong-phase parameters obtained by CLEO. 

PACS numbers: 12.15.Hh, 13.25.Hw, 14.40.Nd 



I. INTRODUCTION 

The angle (f>s (also denoted as 7) is one of the least well- 
constrained parameters of the Unitarity Triangle. The 
measurement that currently dominates sensitivity to 03 
uses B ± — > DK^ decays with the neutral D meson de- 
caying to a three-body final state such as Kgir + n~ [HQ. 
The weak phase appears in the interference between 
b — > cus and b — > ucs transitions. The value of 
03 is determined by exploiting differences between the 
KgTr + ir~ Dalitz plots for D mesons from B + and B~~ 
decay. Theoretical uncertainties in the 03 determination 
in B ± — > DK ± decays are expected to be negligible [H, 
and the main difficulty in its measurement is the very 
low probability of the decays that are involved. How- 
ever, the method based on Dalitz plot analysis requires 
the knowledge of the amplitude of the D° — > Kgir + ir~ 
decay, including its complex phase. The amplitude can 
be obtained from a model that involves isobar and K- 
matrix [4| descriptions of the decay dynamics, and thus 
results in a model uncertainty for the 03 measurement. 
In the latest model-dependent Dalitz plot analyses per- 
formed by BaBar and Belle, this uncertainty ranges from 
3° to 9° @-G3. 

A method to eliminate the model uncertainty using a 



binned Dalitz plot analysis has been proposed by Giri 
et al. Information about the strong phase in the 
D° — > KgTr + ir~ decay is obtained from the decays of 
quantum-correlated D° pairs produced in the ■0(3770) — >• 
D°D° process. As a result, the model uncertainty is re- 
placed by a statistical error related to the precision of the 
strong-phase parameters. This method has been further 
developed in Refs. [ill, EH; where its experimental fea- 
sibility has been shown along with a proposed analysis 
procedure to optimally use the available B decays and 
correlated D° pairs. In this paper, we report the first 
measurement of 03 using a model-independent Dalitz 
plot analysis of the D — > Kgir + Tr~ decay from the mode 
S ± — > DK^, based on a 711 fb _1 data sample (corre- 
sponding to 772 x 10 6 BB pairs) collected by the Belle 
detector at the KEKB asymmetric-energy e + e~ collider. 
This analysis uses the recent measurement of the strong 
phase in D° A'^tt+tt- and D° K^K+K~ decays 
performed by the CLEO collaboration [lj, (3] ■ 



3 



II. THE MODEL-INDEPENDENT DALITZ 
PLOT ANALYSIS TECHNIQUE 



The amplitude of the B+ -> DK+, D 
decay is a superposition of the B + —> D°K 
D°K + amplitudes 



■» K° S TT+1T- 

and B+ - 



A B (m 2 + ,m 2 _) = A + r B e 



i($B+4>3) 



A. 



(1) 



where and m 2 _ are the Dalitz plot variables — the 
squared invariant masses of K%it + and Kgir~ combina- 
tions, respectively, A = A{m\,m 2 _) is the amplitude of 
the TP — > Kgn + 7r~ decay, A = A(m^_, m 2 :) is the ampli- 
tude of the D° — > Kgir + ir~ decay, r B is the ratio of the 
absolute values of the B+ -> D°K+ and B+ -> D°K + 
amplitudes, and S B is the strong-phase difference be- 
tween them. In the case of CP conservation in the D 
decay A(m+,rriL) = A(mt,m 2 h ). The Dalitz plot den- 
sity of the D decay from B + —> DK + is given by 



Pb = \A B \ 2 =\A + r B e 



i(<5s+<fe)^|2 



P + r 2 B P + 2VPP(x + C + y+S), 



(2) 



where P{m 



+ • 1 



L4| 2 , P(m 2 



= \A\ 2 : while 



x + = r B cos(5 B + fa) , y+ = r B sin(S B + fa); (3) 

and the functions C = C(m 2 i _,m 2 _) and S = S(m 2 ^,m 2 _) 
are the cosine and sine of the strong-phase difference 
Soijn 2 ,, m_) = argil — argil between the D Q — > 
K^n + ir~ and D° -> K^tt + it~ amplitudes!]] The equa- 
tions for the charge- conjugate mode B~ —> DK~ are 
obtained with the substitution fa — > —fa and A ^ A; 
the corresponding parameters that depend on the B~ 
decay amplitude are: 

x- = r B cos(S B - </> 3 ) , y- = r B sm(5 B - fa) . (4) 

Using both B charges, one can obtain fa and S B sepa- 
rately. 

Up to this point, the description of the model- 
dependent and model-independent techniques is the 
same. The model-dependent analysis deals directly with 
the Dalitz plot density, and the functions C and S 
are obtained from model assumptions in the fit to the 
D° —> KgTr + ir~ amplitude. In the model-independent 
approach, the Dalitz plot is divided into 27V bins sym- 
metric under the exchange m 2 _ -f-> m\. The bin index 
i ranges from —TV to TV (excluding 0); the exchange 
m 2 + <-> m 2 _ corresponds to the exchange i «-> — i. The 
expected number of events in bin i of the Dalitz plot of 
the D meson from B ± — > DK^ is 



K 



±i 



r B K T* 



2^jK t K-i(x±Ci ±y±Si) 



(5) 



1 This paper follows the convention for strong phases in D decay 
amplitudes introduced in Ref. 1 1 Ijjl . 



where h B is a normalization constant and Ki is the num- 
ber of events in the i th bin of the K { g-K + ir^ Dalitz plot 
of the D meson in a flavor eigenstate. A sample of 
flavor-tagged D° mesons is obtained by reconstructing 
£)*± £) n ± decays (note that charge conjugation is as- 
sumed throughout this paper unless otherwise stated). 
The terms Ci and Sj include information about the func- 
tions C and S averaged over the bin region: 



/ \A\\A\CdV 
j \A\UV $ \A\HV 



(6) 



Here V represents the Dalitz plot phase space and T>i is 
the bin region over which the integration is performed. 
The terms Sj are defined similarly with C substituted by 
S. The absence of CP violation in the D decay implies 
Ci c — i and si s — ^. 

The values of the Cj and Sj terms are measured in the 
quantum correlations of D pairs by charm-factory ex- 
periments operating at the threshold of DD pair pro- 
duction [l3l [l4j]. The measurement involves studies 
of the four-dimensional (4D) density of two correlated 
D —> KgTT + n~ Dalitz plots, as well as decays of a D me- 
son tagged in a CP-eigenstate decaying to K s ir + n~ . The 
wave function of the two mesons is antisymmetric, thus 
the 4D density of two correlated D — > Kgir + ir~ Dalitz 
plots is 



\A cml .( m 2 + ,m 2 _,m' 2 ,m' 2 )\ 2 = \AA' - AA'\ 2 = 



PP 1 + PP' - 2\ / PP'PP'(CC + SS') , 



(7) 



where the primed and unprimed quantities correspond 
to the two decaying D mesons. Similarly, the density of 
the decay D — > KgTr + Tr~, where the D meson is in a 
CP-eigenstate, is 

\A CP (m 2 + ,m 2 _)\ 2 = \A±A\ 2 = P + P ±2> 

CLEO uses these relations to obtain Cj and Sj values. 
Once they are measured, the system of equations ([5]) con- 
tains only three free parameters (x, y, and h B ) for each 
B charge, and can be solved using a maximum likelihood 
method to extract the value of fa. 

We have neglected charm-mixing effects in D decays 
from both the B ± -> DK* process and in quantum- 
correlated DD production. It has been shown [l5| that 
although the charm mixing correction is of first order in 
the mixing parameters xd, yD, it is numerically small (of 
the order 0.2° for xd,Vd ~ 0.01) and can be neglected 
at the current level of precision. Future precision mea- 
surements of fa can account for charm mixing and CP 
violation (both in mixing and decay) using the measure- 
ment of the corresponding parameters. 

In principle, the set of relations defined by Eq. ([5]) can 
be solved without external constraints on Ci and Sj for 
TV > 2. However, due to the small value of r B , there 
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is very little sensitivity to the c, and Sj parameters in 
B ± — > DK^ decays, which results in a reduction in the 
precision on ^3 that can be obtained (Tlj . 

III. CLEO INPUT 

The procedure for a binned Dalitz plot analysis should 
give the correct results for any binning. However, the 
statistical accuracy depends strongly on the amplitude 
behavior across the bins. Large variations of the am- 
plitude within a bin result in loss of coherence in the 
interference term. This effect becomes especially signifi- 
cant with limited statistics when a small number of bins 
must be used to ensure a stable fit. Greater statistical 
precision is obtained for the binning in which the phase 
difference between the D° and D° amplitudes varies as 
little as possible within a bin [l2| . For optimal precision, 
one also has to take the variations of the absolute value 
of the amplitude into account, along with contributions 
from background events. The procedure to optimize the 
binning for the maximal statistical precision of <f>^ has 
been proposed in Ref. ILJ and generalized to the case 
with background in Ref. [lj]. It has been shown that as 
few as 16 bins are enough to reach a statistical precision 
that is only 10-20% worse than in the unbinned case. 

The optimization of binning sensitivity uses the am- 
plitude of the D — > Kgir + Tr~ decay. It should be noted, 
however, that although the choice of binning is model- 
dependent, a poor choice of model results only in a loss 
of precision, not bias, of the measured parameters [l2j]. 
CLEO measured Cj and si parameters for four different 
binnings with M = 8: 

1. Bins equally distributed in the phase difference 
A<5d between the D° and D° decay amplitudes, 
with the amplitude from the BaBar measure- 
ment 0. 

2. Same as option 1. but with the amplitude from the 
Belle analysis |10j |. 

3. Optimized for statistical precision according to the 
procedure from [I2J (see Fig. [T]). The effect of the 
background in B data is not taken into account in 
the optimization. The amplitude is taken from the 
BaBar measurement @. 

4. Same as option 3, but optimized for an analysis 
with high background in B data (e. g., at LHCb). 

Our analysis uses the optimal binning shown in Fig. Q] 
(option 3) as the baseline since it offers better statistical 
accuracy. In addition, we use the equal phase difference 
binning (A^-binning, option 1) as a cross-check. 

The results of the CLEO measurement of Cj and Sj 
for the optimal binning are presented in Table HI The 
same results in graphical form are shown in Fig. [3] The 
values of Cj and calculated from the Belle model [T3| 
are compared to the measurements and are found to be in 




0.5 1 1.5 2 2.5 3 

m 2 (GeV 2 /c 4 ) 



FIG. 1. Optimal binning of the D -»■ ivg7r+7r~ Dalitz plot. 
The color scale indicated corresponds to the absolute value of 
the bin index, \i\. 



TABLE I. Values of a and Si for the optimal binning mea- 
sured by CLEO [3, and calculated from the Belle D — y 
K%ix + tt~ amplitude model. 





CLEO measurement 


Belle model 


Cl 


-0.009 ±0.088 ±0.094 


-0.039 


C2 


±0.900 ±0.106 ±0.082 


±0.771 


C3 


±0.292 ±0.168 ±0.139 


±0.242 


C4 


-0.890 ±0.041 ±0.044 


-0.867 


C5 


-0.208 ±0.085 ±0.080 


-0.246 


C6 


±0.258 ±0.155 ±0.108 


±0.023 


C7 


±0.869 ±0.034 ±0.033 


±0.851 


C8 


±0.798 ±0.070 ±0.047 


±0.662 


Si 


-0.438 ±0.184 ±0.045 


-0.706 


Si 


-0.490 ±0.295 ±0.261 


±0.124 


S3 


-1.243 ±0.341 ±0.123 


-0.687 


SA 


-0.119 ±0.141 ±0.038 


-0.108 


SB 


±0.853 ±0.123 ±0.035 


±0.851 


S 6 


±0.984 ±0.357 ±0.165 


±0.930 


S7 


-0.041 ±0.132 ±0.034 


±0.169 


S8 


-0.107 ±0.240 ±0.080 


-0.596 



reasonable agreement with x 2 = 18.6 for the number of 
degrees of freedom ndf = 16 (the corresponding p-value 
is p = 29%). 

As is apparent from Fig. [21 the chosen binning contains 
bins where the strong phase difference between D° and 
D° amplitudes is close to zero (bins with \i\ = 2, 7, 8) and 
180° (bin with |i| = 4) which provide sensitivity to x±, 
as well as bins with the strong phase difference close to 
90° and 270° (bins with |i| = 1, 3, 5, 6), more sensitive to 
y±. This ensures that the method is sensitive to <f>3 for 
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FIG. 2. Comparison of phase terms (a,Si) for the optimal 
binning measured by CLEO, and calculated from the Belle 
D — > KgTv + -K^ amplitude model. 

any combination of </>3 and 6b values. 



IV. ANALYSIS PROCEDURE 

Equation ([5]) is the key relation used in the analysis, 
but it only holds if there is no background, a uniform 
Dalitz plot acceptance and no cross-feed between bins. 
(Cross-feed is due to invariant-mass resolution and ra- 
diative corrections.) In this section we outline the proce- 
dures that account for these experimental effects. 

A. Efficiency profile 

We note that the Eqs. fl2J, and © do not change 
after the transformation P — > eP when the efficiency pro- 
file e(m+,m?_) is symmetric: e(m+,m?_) = e(m 2 _,m~ i _). 
This implies that if the efficiency profile is the same in all 
of the three modes involved in the measurement (flavor 
D, correlated ^(3770) DD, and D from B -> DK), 
the result will be unbiased even if no efficiency correction 
is applied. 

The effect of non-uniform efficiency over the Dalitz plot 
cancels out when using a flavor-tagged D sample with 
kinematic properties that are similar to the sample for 
the signal B decay. This approach allows for the removal 
of systematic error associated with the possible inaccu- 
racy in the description of the detector acceptance in the 
Monte Carlo (MC) simulation. The center-of-mass (CM) 
D momentum distribution for B — > DK decays is prac- 
tically uniform in the narrow range 2.10 GeV/c < pu < 
2.45 GeV/c. We assume that the efficiency profile de- 
pends mostly on the D momentum and take the flavor- 
tagged sample with an average momentum of po = 2.3 
GeV/c (we use a wider range of D momenta than in 
B DK to increase the statistics). The assumption 



that the efficiency profile depends only on the D momen- 
tum is tested using MC simulation, and the remaining 
difference is treated as a systematic uncertainty. 

While calculating c$ and Sj, CLEO applies an efficiency 
correction, therefore the values reported in their analysis 
correspond to a flat efficiency profile. To use the Cj and 
Si values in the fa analysis, they have to be corrected 
for the Belle efficiency profile. This correction cannot 
be performed in a completely model-independent way, 
since the correction terms include the phase variation 
inside the bin. Fortunately, the calculations using the 
Belle D — s- Kgir + ir~ model show that this correction 
is negligible even for very large non-uniformity of the 
efficiency profile. The difference between the uncorrected 
c» and Si terms and those corrected for the efficiency, 
calculated using the efficiency profile parameterization 
used in the 605 fb _1 analysis [lj|, does not exceed 0.01, 
which is negligible compared to the statistical error. 



B. Momentum resolution 

Momentum resolution leads to migration of events be- 
tween the bins. In the binned approach, this effect can 
be corrected in a non-parametric way. The migration can 
be described by a linear transformation of the number of 
events in bins: 

N'^Ys^kNk, (9) 

where Nk is the number of events that bin k would con- 
tain without the cross-feed, and N[ is the reconstructed 
number of events in bin i. The cross- feed matrix 
is nearly a unit matrix: a.ik <C 1 for i =/= k. It is ob- 
tained from a signal MC simulati on g enerated with the 
amplitude model reported in Ref [lj|. In the case of a 
D —> Kgir + ir~ decay from a B, the cross-feed depends 
on the parameters x and y. However, this is a minor 
correction to an already small effect due to cross-feed; 
therefore it is neglected. 

Migration of events between the bins also occurs due 
to final state radiation (FSR). The Cj and Sj terms in 
the CLEO measurement are not corrected for FSR; we 
therefore do not simulate FSR to obtain the cross-feed 
matrix to minimize the bias due to this effect. Compar- 
ison of the cross-feed with and without FSR shows that 
this effect is negligible. 



C. Fit procedure 

The background contribution has to be accounted for 
in the calculation of the values iV, and Ki. Statistically 
the most effective way of calculating the number of signal 
events (especially in the case of N%, where the statistics is 
a limiting factor) is to perform, in each bin i of the Dalitz 
plot, an unbinncd fit in the variables used to distinguish 
the signal from the background. 
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Two different approaches are used in this analysis to 
obtain the CP violating parameters from the data: sep- 
arate fits in bins, and a combined fit. 

In the first, we fit the data distribution in each bin sep- 
arately, with the number of events for signal and back- 
grounds as free parameters. Once the numbers of events 
in bins N are found, we use them in Eq. [5] to obtain the 
parameters (x± , y± ) . This is accomplished by minimiz- 
ing a negative logarithmic likelihood of the form 

- 2 log C(x, y) = -2j2 logf>((iV;}(x, y),N x , a Ni ), (10) 

i 

where (Ni)(x, y) is the expected number of events in the 
bin i obtained from Eq. [SJ Here, N and are the 
observed number of events in data and the uncertainty 
on Ni, respectively. If the probability density function 
(PDF) p is Gaussian, this procedure is equivalent to a x 2 
fit ; however, the assumption of the Gaussian distribution 
may introduce a bias in the case of low statistics in certain 
bins. 

The procedure described above does not make any as- 
sumptions on the Dalitz distribution of the background 
events, since the fits in each bin are independent. Thus 
there is no associated systematic uncertainty. However, 
in the case of a small number of events and many back- 
ground components this can be a limiting factor. Our 
second approach is to use the combined fit with a com- 
mon likelihood for all bins. The relative numbers of back- 
ground events in bins in such a fit can be constrained 
externally from MC and control samples. In addition, 
for the case of the combined fit, the two-step procedure 
of first extracting the numbers of signal events, and then 
using them to obtain (x, y) is not needed — the expected 
numbers of events (Ni) as functions of (x, y) can be in- 
cluded in the likelihood. Thus the variables (x, y) be- 
come free parameters of the combined likelihood fit, and 
the assumption that the number of signal events has a 
Gaussian distribution is not needed. 

Both approaches are tested with the control samples 
and MC simulation. We choose the combined fit ap- 
proach as the baseline, but the procedure with separate 
fits in bins is also used: it allows a clear demonstration 
of the CP asymmetry in each bin. 

V. EVENT SELECTION 

We use a data sample of 772 x I0 6 BB pairs collected 
by the Belle detector. The decays B ± — > DK^ 1 and 
B ± — >• Dtt^ 1 are selected for the analysis. The neutral 
D meson is reconstructed in the Kgir + ir~ final state in 
all cases. We also select D*^ —> Dn^ decays produced 
via the e + e~ — > cc continuum process as a high-statistics 
sample to determine the Ki parameters related to the 
flavor-tagged D° — > KgTT + ir~ decay. 

The Belle detector is described in detail elsewhere 
[l|| [l7| ■ It is a large-solid-angle magnetic spectrometer 
consisting of a silicon vertex detector (SVD), a 50-layer 



central drift chamber (CDC) for charged particle tracking 
and specific ionization measurement (dE/dx), an array 
of aerogel threshold Cherenkov counters (ACC), time- 
of- flight scintillation counters (TOF), and an array of 
CsI(Tl) crystals for electromagnetic calorimctry (ECL) 
located inside a superconducting solenoid coil that pro- 
vides a 1.5 T magnetic field. An iron flux return located 
outside the coil is instrumented to detect Kl mesons and 
identify muons (KLM). 

Charged tracks are required to satisfy criteria based 
on the quality of the track fit and the distance from 
the interaction point of the beams (IP). We require each 
track to have a transverse momentum greater than 100 
MeV/c, and the impact parameter relative to the IP to 
be less than 2 mm in the transverse and less than 10 mm 
in longitudinal projections. Separation of kaons and pi- 
ons is accomplished by combining the responses of the 
ACC and the TOF with the dE/dx measurement from 
the CDC. Neutral kaons are reconstructed from pairs of 
oppositely charged tracks with an invariant mass M„ n 
within 7 MeV/c 2 of the nominal Kg mass, flight dis- 
tance from the IP in the plane transverse to the beam 
axis greater than 0.1 mm, and the cosine of the angle 
between the projections of Kg flight direction and its 
momentum greater than 0.95. 

The flavor of the neutral D mesons used for Ki deter- 
mination is tagged by the charge of the slow pion in the 
decay D*^ — > D^. The slow pion track is required to 
originate from the D° decay vertex to improve the mo- 
mentum and angular resolution. The selection of signal 
candidates is based on two variables, the invariant mass 

of the neutral D candidates Mq = Mr°-k+-k- an d the 

s i 

difference of the invariant masses of the D*^ and the 
neutral D candidates AM = M^o n + w -^ DW — Al K o 7T + 7T -. 
We retain the events satisfying the following criteria: 
1800 MeV/c 2 < M D < 1920 MeV/c 2 and AM < 150 
MeV/c 2 . We also require the momentum of the D° candi- 
date in the CM frame pu to be in the range 1.8 GeV/c < 
Pd < 2.8 GeV/c to reduce the effect of the efficiency pro- 
file on the (f>3 measurement (see Sec. IIV A]) . About 15% 
of selected events contain more than one -D candidate 
that satisfies the requirements above; in this case we keep 
only one randomly selected candidate. 

Selection of B ± — > DK^ and B ± -> Z?7r ± samples 
is based on the CM-cnergy difference AE = Ei — 
£"bcam and the beam-constrained B meson mass M^ c = 
V^bcam "(Ek) 2 : where E hcam is the CM beam energy, 
and Ei and pi are the CM energies and momenta of the B 
candidate decay products. We select events with Mb c > 
5.2 GeV/c 2 and \AE\ < 0.18 GeV for further analysis. 
We also impose a requirement on the invariant mass of 
the neutral D candidate \M K o v + v - —M D o\ < 11 MeV/c 2 . 

Further separation of the background from e + e~ — > qq 
(q = u, d, s, c) continuum events is done by calculating 
two variables that characterize the event shape. One is 
the cosine of the thrust angle cos#thri where # t hr is the 
angle between the thrust axis of the B candidate daugh- 
ters and that of the rest of the event, calculated in the 
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CM frame. The other is a Fisher discriminant T com- 
posed of 11 parameters [llj]: the production angle of the 
B candidate, the angle of the B thrust axis relative to 
the beam axis, and nine parameters representing the mo- 
mentum flow in the event relative to the B thrust axis 
in the CM frame. We use the AE, Mb c , cos6*thr, and T 
variables in the maximum likelihood fit. 

In both flavor D° and B ± -> DK* (B ± -> L>tt ± ) 
samples, the momenta of the tracks forming a D candi- 
date are constrained to give the nominal D° mass in the 
calculation of the Dalitz plot variables. 



VI. FLAVOR-TAGGED SAMPLE D ->■ , 

D -> K%-K + it- 

The number of events Ki in bin i of the flavor-tagged 
D — > KgTT + ir~ decay is obtained from a two-dimensional 
unbinned fit to the distribution of Mo and AM variables. 
The fits in each Dalitz plot bin are performed indepen- 
dently. The fit uses a signal PDF and two background 
components: purely random combinatorial background 
and background with a real D° and random slow pion 
track. The signal distribution is a product of the PDFs 
for Md (triple Gaussian) and AM (sum of bifurcated 
Student's i-distribution and bifurcated Gaussian distri- 
bution) . The combinatorial background is parameterized 
by a linear function in Mo and by a function with a kine- 
matic threshold at the ir + mass in A M: 

Pcomb(AM) = y/y(l + Ay[l + B{M D - m D o)}) e - yC 

(11) 

where y = AM — m v + , and m^o are the nominal 
masses of n + and D°, respectively, and A, B, and C are 
free parameters. A small correlation between the Mo 
and AM distributions is introduced that is controlled by 
the parameter B. The random slow pion background is 
parameterized as a product of the signal Mo distribution 
and combinatorial AM background shape. 

The parameters of the signal and background distribu- 
tions are obtained from the fit to data. The parameters 
of the signal PDF are constrained to be the same in all 
bins. The free parameters in each bin are the number of 
signal events Ki, the parameters of the background dis- 
tribution, and fractions of the background components. 

The fit results from the flavor-tagged D sample inte- 
grated over the whole Dalitz plot are shown in Fig. [3] 
The number of signal events calculated from the integral 
of the signal distribution is 426938 ±825, the background 
fraction in the signal region \Mjj — m D o\ < 11 MeV/c 2 , 
144.5 MeV/c 2 < AM < 146.5 MeV/c 2 is 10.1 ± 0.1%. 
The signal yield in bins is shown in Table |TTJ 




H D (GeV/c ! ) AM (GeV/c 2 ) 

FIG. 3. Projections of the flavor-tagged D — > Dir^ , 
D -> A-gvr+Tr- data with 1.8 GeV/c < p D < 2.8 GeV/c. (a) 
M D distribution for 144.5 MeV/c 2 < AM < 146.5 MeV/c 2 . 
(b) AM distribution for 1854 MeV/c 2 < M D < 1876 MeV/c 2 . 
Histograms show the fitted signal and background contri- 
butions, points with the error bars are the data. The full 
D — > K%-k + -k~ Dalitz plot is used. 



TABLE II. Signal yields in Dalitz plot bins for the flavor- 
tagged D* ± — > Dir , D — > K%-k + h~ sample with 
1.8 GeV/c <p D < 2.8 GeV/c. 



Bin i 


Ki 


K-i 


1 


43261 ± 255 


8770 ±124 


2 


58005 ± 268 


1827 ± 63 


3 


62808 ± 274 


1601 ± 58 


4 


44513 ± 253 


26482 ± 202 


5 


21886 ± 177 


13146 ± 143 


6 


28876 ± 197 


1765 ± 68 


7 


48001 ± 265 


22476 ± 196 


8 


9279 ± 125 


26450 ± 181 


Total 


426938 ± 825 



VII. SELECTION OF B ± -> Dvr* AND B ± -> DK ± 
SAMPLES 

The decays B^ DK^ and B ± — ^ Z?7r ± have simi- 
lar topology and background sources and their selection 
is performed in a similar way. The mode B ± -> £)tt ± 
has an order of magnitude larger branching ratio and a 
smaller amplitude ratio rs ~ 0.01 due to the ratio of 
weak coefficients \V ub V* d \/\V cb V* d \ - 0.02 and the color 
suppression factor. This results in the small CP viola- 
tion in this mode, therefore it is used as a control sample 
to test the procedures of the background extraction and 
Dalitz plot fit. In addition, signal resolutions in AE and 
A/bc and the Dalitz plot structure of some background 
components are constrained from the control sample and 
used in the signal fit. 

The number of signal events is obtained by fitting the 
4D distribution of variables Mb c , AE, cos6* t hr and T. 
The fits to the B ± — > Z?7r ± and B ± — > DK^ samples use 
the following three background components in addition 
to the signal PDF: 

• Combinatorial background from the process 
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e + e — > qq, where q — (u, d, s, c). 

• Random BB background, in which the tracks form- 
ing the i? ± — > Dtt^ 1 candidate come from decays 
of both B mesons in the event. The number of 
possible B decay combinations that contribute to 
this background is large, therefore both the Dalitz 
distribution and (Mh c ,AE) distribution are quite 
smooth. 

• Peaking BB background, in which all tracks form- 
ing the B^ — > candidate come from the same 
B meson. This kind of background is dominated 
by B — > D*ir decays reconstructed without the tt 
or 7 from the D* decay. 

In addition, the B^ — > DK^ fit includes a fourth com- 
ponent that models B ± — > Dn^ decays in which the pion 
is misidentified as a kaon. 

Thc PDF for the signal parameterization (as well as for 
each of the background components) is a product of the 
(M bc ,AE) and (cos thr , F) PDFs. Thc (M hc , AE) PDF 
is a 2D double- Gaussian function, which has a correlation 
between M^ c and AE. The double- Gaussian function 
models both the core and tails of the distribution. The 
(cos^thr,-?") distribution is parameterized by the sum of 
two functions (with different coefficients) of the form 

p(x,F) =exp(Cia; + C 2 x 2 + C 3 x 3 )x 

G(F, F (x), <t fl (x), a FR (x)), 

where x = cos# t hr, G(F, F,aL,an) is thc bifurcated 
Gaussian distribution with the mean F and the widths 
<7l and <tr, and functions Fq, op^ and <jfr are polynomi- 
als that contain only even powers of x. The parameters 
of the signal PDF are obtained from the signal MC simu- 
lation. However, to account for thc possible imperfection 
of the simulation, we allow all the width parameters to 
scale by a common factor, which is obtained from the 
B ± — > Dir^ sample. 

The combinatorial background from continuum 
e + e~ — > qq production is obtained from the experimen- 
tal sample collected at a CM energy below the T(4S*) 
resonance (off- resonance data). The parameterization in 
variables (cos £> t h r , J 7 ) follows Eq. (|T2"j) . The parameter- 
ization in (Afbc, AE) is the product of an exponential 
distribution in AE and the empirical shape proposed by 
the ARGUS collaboration [H] in M hc : 

Pcomb(M hc , AE) = exp(-aAE)M hcy /yexp(-cy), (13) 

where y = l-M hc /E hcam , E hc&m is the CM beam energy, 
and a and c are empirical parameters. 

The parameters for random and peaking BB back- 
grounds are obtained from a generic MC sample. Gener- 
ator information is used to distinguish between the two: 
the latter contains only the events in which the candidate 
is formed of tracks coming from both B mesons. The 
(Mb c ,A£) distributions for each of these backgrounds 
are parameterized by the sum of three components: 



• the product of an exponential (in AE) and Argus 
(in Mbc) functions, as for continuum background 
(as expected, this component dominates the ran- 
dom BB background); 

• the product of an exponential in the AE and bi- 
furcated Gaussian distribution in Mb c , where the 
mean of the Gaussian distribution is linear as a 
function of AE; and 

• a two-dimensional Gaussian distribution in AE and 
Mbc, which includes a correlation and is asymmet- 
ric in Mhc- This component is small compared to 
the random BB contribution, but dominates the 
peaking BB background, which mostly consists of 
partially reconstructed B decays. 

The peaking background coming from B + B~ and 
B°B° decays is treated separately in (M\> C ,AE) vari- 
ables, while a common (cos 0thr, F) distribution is used. 
In the case of the B ± -> fit, B ± -> events 

with the pion misidentified as a kaon are treated as 
a separate background category. The distributions of 
A/bc,A_B and co&Ot^^F variables are parameterized in 
the same way as for the signal events and are obtained 
from MC simulation. 

The Dalitz plot distributions of the background com- 
ponents are discussed in the next section. Note that the 
Dalitz distribution is described by the relative number of 
events in each bin. The numbers of events in bins can 
be free parameters in the fit, thus there will be no uncer- 
tainty due to the modeling of the background distribution 
over the Dalitz plot in such an approach. This procedure 
is justified for background that is either well separated 
from the signal (such as peaking BB background in the 
case of B^ — > Dir^), or is constrained by a much larger 
number of events than the signal (such as the continuum 
background) . 

The results of thc fit to B ± and B ± 

DK^ data with the full Dalitz plot taken are shown 
in Figs. @] and O respectively. We obtain a total of 
19106 ± 147 signal B ± Dn ± events and 1176 ± 43 
signal B ± — > DK^ events — 55% more than in the 605 
fb -1 model-dependent analysis [lo| . The improvement 
partially comes from the larger integrated luminosity of 
the sample, and partially from the larger selection effi- 
ciency due to improved track reconstruction. 



VIII. DATA FITS IN BINS 

The data fits in bins for both B ± — s- £>7r ± and B ± -> 
DK samples are performed with two different proce- 
dures: separate fits for the number of events in bins and 
the combined fit with the free parameters (x, y) as dis- 
cussed in Sec. IIVC1 The combined fit is used to obtain 
the final values for {x,y) 1 while the separate fits provide 
a crosscheck of the fit procedure and a way to visualize 
the extent of CP violation within the sample. A study 
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FIG. 4. Projections of the B ± Dvr* data, (a) M bc dis- 
tribution with \AE\ < 30 MeV and cos^thr < 0.8 require- 
ments, (b) AE distribution with M hc > 5270 MeV/c 2 and 
cos^thr < 0.8 requirements, (c) cos# t h r and (d) T distribu- 
tions with \AE\ < 30 MeV and M bc > 5270 MeV/c 2 require- 
ments. Histograms show the fitted signal and background 
contributions, points with error bars are the data. The entire 
D — > KgTr + Tr~ Dalitz plot is used. 
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FIG. 5. Projections of the B ± -s> DK^ data, (a) M hc dis- 
tribution with \AE\ < 30 MeV and cos^thr < 0.8 require- 
ments, (b) AE distribution with M hc > 5270 MeV/c 2 and 
cos#thr < 0.8 requirements, (c) cos# t hr and (d) T distribu- 
tions with \ AE\ < 30 MeV and M hc > 5270 MeV/c 2 require- 
ments. Histograms show the fitted signal and background 
contributions, points with error bars are the data. The entire 
D Kgir + ir~ Dalitz plot is used. 



with MC pseudo-experiments is performed to check that 
the observed difference in the fit results between the two 
approaches agrees with expectation. 

In the case of separate fits in bins, we first perform 
the fit to all events in the Dalitz plot. The fit uses back- 
ground shapes fixed to those obtained from fits to the 
generic MC samples of continuum and BB decays. The 
signal shape parameters are fixed to those obtained from 
a fit to the signal MC sample except for the mean value 
and width scale factors of AE and Mb c PDFs. As a next 
step, we fit the 4D (Mb c , AE, cos (?thr, J~) distributions in 
each bin separately, with the signal peak positions and 
width scale factors fixed to the values obtained from the 
fit to all events. The free parameters of each fit are the 
number of signal events, and the number of events in each 
background category. 

The numbers of signal events in bins for the B ± — > 
Z?7r ± sample extracted from the fits are given in Table Hill 
These numbers are used in the fit to extract (x, y) using 
Eq. ([5]) after the cross-feed and efficiency correction for 
both Ni and Ki . Figure [5] illustrates the results of this 
fit. The numbers of signal events in each bin for B + and 
B~ are shown in Fig. Eta) together with the numbers 
of events in the flavor-tagged D° sample (appropriately 
scaled). The difference in the number of signal events 
shown in Fig. Etb) does not reveal CP violation. Fig- 
ures [51c) [(d)] show the difference between the numbers 
of signal events for B + [B~] data and scaled flavor-tagged 
D° sample, both for the data and after the (x, y) fit. The 
X 2 /ndf is reasonable for both the (x,y) fit and the com- 
parison with the flavor-specific CP conserving amplitude. 

Unlike B ± — s- Z?7r ± , the B ± -t DK^ sample has sig- 
nificantly different signal yields in bins of B + and B~ 
data (see Fig. \T[b) and Table HVJ). The probability to 
obtain this difference as a result of a statistical fluctu- 
ation is 0.42%. This value can be taken as the model- 
independent measure of the CP violation significance. 
The significance of 03 being nonzero is in general smaller 
since (j>3 7^ results in a specific pattern of charge asym- 
metry. The fit of the signal yields to the expected pat- 
tern described by the parameters (x, y) is of good qual- 
ity [Figs. El(c,d)], which is consistent with the hypothesis 
that the observed CP violation is solely explained by the 
mechanism involving nonzero ())■$. 

The default combined fit constrains the random BB 
background in bins from the generic MC, and takes the 
(x±,y±) variables as free parameters. Fits to B + and 
B~ data are performed separately. The plots illustrat- 
ing the combined fit results are given in the Appendix. 
Additional free parameters are the yields of continuum 
and peaking BB backgrounds in each bin, the fraction 
of the random BB background, and the means and scale 
factors of the signal M^ c and AE distributions. The val- 
ues of (x, y) are then corrected for the fit bias obtained 
from MC pseudo-experiments. The value of the bias de- 
pends on the initial x and y values and is of the order 
5 x 10~ 3 for the B^ — > DK ± sample and less than 1CT 3 
for B^ — > Dtt^ sample. 
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TABLE III. Signal yields in Dalitz plot bins for the B ± -> TABLE IV. Signal yields in Dalitz plot bins for the B ± -> 
D% ± , D —> K% Ti + n sample with the optimal binning. DK , D — > Jfg7r + 7r sample with the optimal binning. 



Bin i 




N? 




Bin i 






-8 


564.2 ± 25.3 


587.0 ± 25.7 




-8 


49.8 ± 8.2 


37.8 ± 7.5 


-7 


462.3 ± 23.8 


462.8 ± 23.9 




-7 


42.2 ± 8.6 


24.9 ± 7.2 


-6 


47.9 ± 7.7 


39.2 ± 7.2 




-6 


0.0 ± 1.9 


3.4 ±2.9 


-5 


314.1 ± 19.0 


286.2 ± 18.2 




-5 


9.6 ±4.5 


23.6 ± 6.2 


-4 


592.6 ± 26.5 


645.7 ± 27.8 




-4 


32.9 ±7.5 


42.1 ± 8.3 


-3 


22.2 ±6.2 


27.2 ±6.3 




-3 


3.5 ±2.8 


0.7 ±2.5 


-2 


42.7 ±7.6 


54.0 ±8.7 




-2 


11.3 ±4.1 


0.0 ± 1.3 


-1 


190.8 ± 15.4 


210.8 ± 16.3 




-1 


16.6 ± 5.4 


7.7 ±4.4 


1 


959.2 ± 32.6 


980.2 ±33.1 




1 


37.6 ±8.0 


65.1 ±9.9 


2 


1288.7 ± 37.0 


1295.9 ±37.1 




2 


68.6 ±9.6 


75.5 ±9.8 


3 


1395.8 ± 38.4 


1352.2 ±37.9 




3 


83.4 ± 10.1 


82.4 ± 10.2 


4 


1045.5 ± 34.7 


1065.1 ±34.9 




4 


49.3 ± 9.1 


86.5 ± 11.4 


5 


479.3 ± 23.3 


532.2 ±24.5 




5 


34.0 ± 7.3 


38.3 ± 7.6 


6 


623.7 ±26.0 


663.5 ±26.7 




6 


34.8 ± 6.8 


41.9 ±7.5 


7 


1081.0 ±35.3 


1049.2 ± 34.8 




7 


70.8 ± 10.6 


46.4 ± 9.0 


8 


210.0 ± 16.1 


212.1 ± 16.3 




8 


9.4 ±4.3 


14.2 ± 5.1 


Total 


9467.1 ± 103.6 


9639.1 ± 104.7 




Total 


574.9 ± 29.9 


601.6 ±30.8 
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FIG. 6. Results of the fit to the B -> Dtt* sample, (a) Sig- 
nal yield in bins of the D — > K%tt + tt~ Dalitz plot: from 
B~ — > Dir~ (red), B + — > Dtv + (blue) and flavor sam- 
ple (histogram), (b) Difference of signal yields between the 
B + — > Dn + and B~ — > Dtv~ decays, (c) Difference of signal 
yields between the B~ — > Dn~ and flavor samples (normal- 
ized to the total B~ — > Dn~ yield): yield from the separate 
fits (points with error bars), and as a result of the combined 
(x,y) fit (horizontal bars), (d) Same as (c) for B + — >• Dn + 
data. 



FIG. 7. Results of the fit of the B DK* sample, (a) Sig- 
nal yield in bins of the D — > Kgir + ir~ Dalitz plot: from 
B~ -> DK~ (red), B + -> DK + (blue) and flavor sam- 
ple (histogram), (b) Difference of signal yields between the 
B + — > DK + and B~ — > DK~ decays, (c) Difference of signal 
yields between the B~ — > DK~ and flavor samples (normal- 
ized to the total B~ —¥ DK~ yield): yield from the separate 
fits (points with error bars), and as a result of the combined 
(x,y) fit (horizontal bars), (d) Same as (c) for B + —¥ DK + 
data. 
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TABLE V. (x,y) parameters and their statistical correlations from the combined fit of the B — > Dir ± and B — > DK ± 
samples. The first error is statistical, the second error is systematic, and the third error is due to the uncertainty on the d and 
Si parameters. 



Parameter 




— > DK ± 


X- 

y- 

corr(x_,y_) 


-0.0045 ± 0.0087 ± 0.0049 ± 0.0026 
-0.0231 ± 0.0107 ± 0.0041 ± 0.0065 
-0.189 


+0.095 ± 0.045 ± 0.014 ± 0.010 
+0.137lg;^ ± 0.015 ± 0.023 
-0.315 


x+ 

y+ 

corr(x+,y + ) 


-0.0172 ± 0.0089 ± 0.0060 ± 0.0026 
+0.0129 ± 0.0103 ± 0.0059 ± 0.0065 
-0.205 


-0.110 ± 0.043 ± 0.014 ± 0.007 
-0.050l2 o55 ± 0.011 ± 0.017 
+0.059 
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FIG. 8. One, two, and three standard deviation (x, y) confi- 
dence levels for (a) B — > Dn ± and (b) B — > DK ± decays 
(statistical only). Note different scale for B — > D-k^ and 
B ± — > DK ± modes. The weak phase 03 appears as half the 
opening angle between (x + ,y + ) and (x-,y-) vectors. 



The values of (x, y) parameters and their statistical 
correlations obtained from the combined fit for B^ — s- 
Dir^ control sample and B^ — > DK^ sample are given in 
Table El The measured values of {x±,y±) for both sam- 



ples with their statistical likelihood contours are shown 
in Fig. E 



IX. SYSTEMATIC ERRORS 

Systematic errors in the (x, y) fit are obtained for the 
default procedure of the combined fit with the optimal 
binning. The systematic errors are summarized in Ta- 
ble EH 

The uncertainty due to the signal shape used in the fit 
includes the sources listed below: 

1. The choice of parameterization used to model 
the shape. The corresponding uncertainty is es- 
timated by using the non-parametric Keys PDF 
function [201 ] instead of the parameterized distri- 
bution. 

2. Correlation between the (Mt, c , AE) and 
(cos #thr, J~) distributions. To estimate its ef- 
fect, we use a 4D binned histogram to describe the 
distribution. 

3. The MC description of the (cos#thr, -T 7 ) distribu- 
tion. Its effect is estimated by floating the parame- 
ters of the distribution in the fit to the B ± — > Dw^ 
control sample. 

4. The dependence of the signal width on the Dalitz 
plot bin. The uncertainty due to this effect is esti- 
mated by performing the B ± — > Z?7r ± fit with the 
shape parameters floated separately for each bin, 
and then using the results in the fit to B ± — > DK^ 
data. 

We do not assign an uncertainty due to the difference 
in (Mbc, AE) shape between the MC and data since the 
width of the signal distribution is calibrated using B ± — > 
Z?7r ± data. 

For the uncertainty due to the continuum background 
shape, we use the same four sources as considered for the 
signal distribution. The uncertainty due to the choice 
of parameterization is estimated by using the Keys PDF 
as an alternative. The effect of correlation between the 
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(Mbc, AE) and (cos 0thr, J~) distributions is estimated by 
using a distribution split into the sum of two compo- 
nents (u, <i, s and charm contributions) with independent 
(Mbc, AE) and (cos t hr, J 7 ) shapes. The uncertainty due 
to the MC description of the (Mbc AE) and (cos 6 t hr, J 7 ) 
distributions is estimated by floating their parameters in 
the B^ — > Dtt^ 1 fit. To estimate the effect of correla- 
tion of the shape with the Dalitz plot variables we fit the 
shapes separately in each Dalitz plot bin. 

The uncertainties due to the shapes of random and 
peaking BB backgrounds arc estimated differently for 
the B ± -> Dn ± and B ± -> DK± samples. In the B ± -> 
Dtt^ case, the effect of the background shape uncertainty 
is estimated by performing a fit with AE > —0.1 GeV: 
this requirement rejects the peaking BB background and 
a large part of the random BB background. In the case 
of the B ± — > DK^ sample, the uncertainty is estimated 
by performing an alternative fit with the (Md C , AE) and 
(cos #thrj F) shapes taken from the B ± — > Dn^ sample 
(shifted by 50 McV in AE to account for difference in 
pion and kaon masses in the calculation of AE) . 

In the case of the fit to the B ± — > DK^ sample, the 
uncertainty due to the B — > Dtt^ background shape 
in the (cos #thr, F) variables is estimated by taking the 
(cos#thr, J 7 ) shape for signal events. The Dalitz plot dis- 
tribution uncertainty is estimated by using the number 
of flavor-tagged events in bins, rather than the number of 
£? ± — > Dn^ events used in the default fit. Uncertainties 
due to correlations are treated in the same way as in the 
case of the signal distribution. 

There is an uncertainty due to the Dalitz plot efficiency 
shape because of the difference in average efficiency over 
each bin for the flavor and B^ — > DK^ samples. A 
maximum difference of 1.5% is obtained in a MC study. 
The uncertainty is taken to be the maximum of the two 
comparable in size quantities: 

• the root mean square (RMS) of x and y from smear- 
ing the numbers of events in the flavor sample Ki 
by 1.5% (larger for y parameters); 

• the bias of x and y between fits with and without ef- 
ficiency correction for Ki obtained from signal MC 
(larger for x parameters). 

The uncertainty due to cross-feed of events between 
bins is estimated by varying the momentum resolution by 
20% — the MC resolution scaling factor obtained from 
the fit to i? ± — > Dn^ sample — and by taking the bias 
between the fits with and without final state interactions 
taken into account. 

The uncertainty arising from the finite sample of flavor- 
tagged D —> Kgir + ir~ decays is evaluated by varying the 
values of Ki within their statistical errors. 

The final results for (x, y) are corrected for the fit bias 
obtained from fits of MC pseudo-experiments. The un- 
certainty due to the fit bias is taken from the difference 
of biases for various input values of x and y. 

The uncertainty due to the limited precision of the Ci 
and Si parameters is obtained by smearing the Ci and Sj 



values within their total errors and repeating the fits for 
the same experimental data. We have performed a study 
of this procedure using both MC pseudo-experiments and 
analytical calculations. We find that the uncertainty ob- 
tained in this way is sample-dependent for small B data 
samples and its average value scales in inverse proportion 
to the square root of the sample size. It reaches a con- 
stant value for large B data samples (in the systematics- 
dominatcd case). This explains the somewhat higher un- 
certainty compared to the CLEO estimate given in flij ]. 
which was obtained in the limit of a very large B sample. 
In addition, the uncertainty in (x, y) is proportional to 
r\e, and, thus, the uncertainty in the phases fa and 6 b is 
independent of Tb- As a result, the uncertainty of (x,y) 
in the B^ — > DK^ sample fit is 3-4 times larger than in 
the — > Dtt^ 1 sample. 



X. RESULTS FOR fa, r B AND S B 

We use the fre quen tist approach with the Feldman- 
Cousins ordering [2l| to obtain the physical parame- 
ters fi = (fa^rs^s) from the measured parameters 
z = (x , y '_, x+, y+), as was done in previous Belle anal- 
yses H,[l3l. In essence, the confidence level a for a set of 
physical parameters /i is calculated as 



a{n) 



p(z\fj,)dz J I p(z\n)dz , 



(14) 



where p(z\fi) is the probability density to obtain the mea- 
surement result z given the set of physics parameters /z. 
The integration domain T>(fi) is given by the likelihood 
ratio (Feldman-Cousins) ordering: 



p{z\ii) 



> 



p(zo\n) 



p(z\fJ,best(z)) Pizolllbestizo))' 



(15) 



where Hi,est(z) is [i that maximizes p(z\fx) for the given 
z, and zq is the result of the data fit. 

In contrast to previous Belle analyses [l,[l(|, the proba- 
bility density p(z\fx) is a multivariate Gaussian PDF with 
the errors and correlations between x± and y± taken from 
the data fit result. In the previous analyses, this PDF was 
taken from MC pseudo-experiments. 

As a result of this procedure, we obtain the confidence 
levels (CL) for the set of physical parameters ^3, t\b, and 
8b ■ The confidence levels for one and two standard devia- 
tions are taken at 20% and 74% (appropriate for the case 
of a three-dimensional Gaussian distribution) . The pro- 
jections of the 3D surfaces bounding one, two, and three 
standard deviations volumes onto (^3,^) and (fa, 6b) 
planes are shown in Fig. [5] 

Systematic errors in [i are obtained by varying the mea- 
sured parameters z within their systematic errors (Gaus- 
sian distributions are assumed) and calculating the RMS 
of /Zbest(z)- In this calculation we assume that the sys- 
tematic errors arc uncorrclatcd between the B + and B~ 
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TABLE VI. Systematic uncertainties in the (x, y) measurements for B — > Dir ± and B — > DK ± samples in units of 10 







B ± — 


► Dtt* 






DK^ 1 




Source of uncertainty 


Ax- 


Ay- 


Ax+ 


Ay+ 


Ax- 


Ay. 


Ax+ 


Ay+ 


Signal shape 


0.9 


1.9 


1.1 


5.0 


7.3 


7.4 


7.3 


5.1 


u, d, s, c continuum background 


0.9 


1.4 


0.8 


1.3 


6.7 


5.6 


6.6 


3.2 


BB background 


3.3 


1.6 


4.5 


1.1 


7.7 


8.4 


7.4 


5.4 


B — > Dir^ background 


— 


— 


— 


— 


1.2 


4.2 


1.9 


1.9 


Dalitz plot efficiency 


3.0 


1.9 


3.2 


1.6 


4.8 


2.0 


5.6 


2.1 


Crossfeed between bins 


0.3 


0.6 


0.1 


0.7 


0.0 


3.9 


0.1 


1.0 


Flavor-tagged statistics 


1.7 


2.0 


1.6 


2.0 


1.5 


2.7 


1.7 


1.9 


Fit bias 


0.4 


0.5 


0.4 


0.5 


3.2 


5.8 


3.2 


5.8 


d and Si precision 


2.6 


6.5 


2.6 


6.5 


10.1 


22.5 


7.2 


17.4 


Total without Ci,Si precision 


4.9 


4.1 


6.0 


5.9 


14.0 


15.3 


14.1 


10.6 


Total 


5.6 


7.7 


6.5 


8.8 


17.3 


27.2 


15.9 


20.4 




number of MC pscudo-expcrimcnts. The statistical sig- 
nificance equals 99.64% or 2.9 standard deviations. This 
value is in good agreement with the x 2 probability from 
the difference of the number of events in bins for B + and 
B~ data. With the systematic uncertainties included, 
the significance decreases to 99.35% or 2.7 standard de- 
viations. 



XI. CONCLUSION 



FIG. 9. Two-dimensional projections of confidence region 
onto {4>3,Sb) and (03, re) planes (one, two, and three stan- 
dard deviations, statistical only). 



samples. In the case of Cj and Si systematics, we test this 
assumption. When the fluctuation in Cj and Sj is gener- 
ated, we perform fits to both B + and B~ data with the 
same fluctuated (ci, s{). We observe no significant corre- 
lation between the resulting X- and x + (y_ and y+). 
The final results are: 

3 = (77.3± 14 ^±4.1±4.3)° 

r B =0.145 ±0.030 ±0.010 ±0.011 (16) 

8 B = (129.9 ± 15.0 ± 3.8 ± 4.7)°, 

where the first error is statistical, the second is systematic 
error without Ci and Si uncertainty, and the third error is 
due to Ci and Sj uncertainty. Extraction of 03, Tb, and 8 b 
has a two-fold ambiguity, (0 3 , 8) and (0 3 + 18O°, (5+180°), 
leading to the same values of x± and y± . Here we choose 
the solution that satisfies < 03 < 180°. 

The significance of CP violation (03 being non-zero) 
is calculated as the CL of the point 03 = 0. This calcula- 
tion accounts for a small deviation from Gaussian errors 
for x and y observed and parameterized using a large 



We report the results of a measurement of the Unitar- 
ity Triangle angle 03 using a model-independent Dalitz 
plot analysis of D — > Kgir + ir~ decay in the process 
B — > DK . The measurement was performed with 
the full data sample of 711 fb" 1 (772 x 10 6 BB pairs) 
collected by the Belle detector at the T(4S) resonance. 
Model independence is achieved by binning the Dalitz 
plot of the D — > KgTT + TT~ decay and using the strong- 
phase coefficients for bins measured by the CLEO ex- 
periment [3]. We obtain the value 3 = (77.3j^4;g ± 
4.1 ± 4.3)°; of the two possible solutions we choose the 
one with < 03 < 180°. We also obtain the value of 
the amplitude ratio r B = 0.145 ± 0.030 ± 0.010 ± 0.011. 
In both results, the first error is statistical, the second is 
systematic error without Cj and Si uncertainty, and the 
third error is due to Cj and Sj uncertainty. 

This analysis is the first application of a novel 
method for measuring 03. Compared to the result of 
the model-dependent measurement performed by Belle 
with the B ± -> L»if ± mode, 3 = (80.8±i|"| ± 
5.0(syst)±8.9(model))° [l|, this measurement has some- 
what poorer statistical precision despite a larger data 
sample used. There are two factors responsible for lower 
statistical sensitivity: 1) the statistical error for the same 
statistics is inversely proportional to the tb value, and 
the central value of tb in this analysis is smaller, and 
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2) the binned approach is expected to have the statisti- 
cal precision that is, on average, 10-20% poorer than the 
unbinned one 

More important is that the large model uncertainty 
of the model-dependent result (8.9°) is replaced by the 
purely statistical uncertainty of 4.3° due to the limited 
size of the CLEO -0(3770) data sample. Although the 
model-independent approach does not offer significant 
improvement over the unbinned model-dependent Dalitz 
plot analysis with the current data sample, it is promi sing 
for future measurements at super flavor factories [22l |23| 
and LHCb (24[. We expect that the statistical error of 
the </>3 measurement using the statistics of a 50 ab^ 1 
data sample that will be available at a super-B factory 
will reach 1-2°. With the use of BES-III data [H the 
error due to the phase terms in the D — > Kgir + ir~ decay 
will decrease to 1° or less. We also expect that the exper- 
imental systematic error can be kept at the level below 
1°, since most of its sources arc limited by the statistics 
of the control channels. 
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APPENDIX 

The results of the combined fit to B + — > DK + and 
B~ — > DK~ samples separately for each bin of the Dalitz 
plot are shown in Figs. ITUl and [TTl respectively. The plots 
show the projections of the data and the fitting model 
on the AE variable, with the additional requirements 
M hc > 5270 MeV/c 2 and cos0 th r < 0.8. 
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FIG. 11. Projections of the combined fit of the B — > DK sample on AE for each Dailtz plot bin, with the M bc > 5270 MeV/c 2 
and cos # t h r < 0.8 requirements. The fill styles for the signal and background components are the same as in Fig. \S\ 



